scipy optimize.root_scalar() with multiple arguments each being a list